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Abstract 

We present a model for neutrino flares in accreting black holes based on the 
injection of a non-thermal population of relativistic particles in a magnetized corona. 
The most important products of hadronic and photohadronic interactions at high 
energies are pions. Charged pions decay into muons and neutrinos; muons also decay 
yielding neutrinos. Taking into account these effects, coupled transport equations 
are solved for all species of particles and the neutrino production is estimated for 
the case of accreting galactic black holes. 
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1 Introduction 



The existence of a broad-band X-ray/soft gamma-ray spectrum of accret- 
ing black holes in binary systems strongly suggests the presence of a very 
hot plasma (corona) around the central object. This plas ma Comptonizes the 



soft X -ray photons produced in the inner accretion disk (jShakura &: Sunyaev 



19731 ) generating a power-law, high-energy feature in the spectrum. A wide 



variety of stable two-temperature corona models have been developed since 
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Shapiro et al.l (119761 ) introduced a temperature differenciation for ions and 
electrons in the inner acc retion flow of Cygnus X-1 (see, for instance, the 
review by iPoutanen 1 119981 ) . 



Galactic black holes are found in different spectral states. The most typical 
are the high-soft, in which the spectrum is dominated by thermal emission 
of the disk, and the low-hard, characterized by the presence of a hot corona 
around the compact object and steady jets. The existence of these jets is 
supported by observational e vidence of radio emission in almost every black 
hole binary in the hard state fICorbel et aDl2000l : IPender fc Bellonill2004h . The 
radio emission is generally too bright to be produced by thermal electrons in 
the accretion flow. 



A strong correl ation has been found between the radio and X-ray luminosity 
in the hard state ( Corbel et al. 20031 ). This suggests that the X-ray emission 
is strongly coupled to that of the jets (e.g., iFalcke et al.ll2004l ). Since the jets 
flow out of the corona and are thus tightly associated to it, a model in which 
the X-ray emission is mostly from the disk -l-corona and radio is fr om the 
jets is compatible with the observations (INaravan fc McClintockll2008[) . Other 
possibilities are discussed, for instance, by iMarkoff et al.l (120031 ). 



In two-temperature magnetized plasmas around black holes it is reason- 
able to expect significant deviations from purely Maxwellian distributions for 
the particles. Such deviations are the result of the injectio n of non-thermal 
populations of particles. Recently, IVurm &: PoutanenI (120091 ) have presented a 
model for time- dependent non-therna al (purely leptonic) emission from mag- 
netized coronae. Romero et al. (l2010h have studied the case of a steady state 
magnetized corona with injection of both relativistic electrons and protons. 
In the present work we extend the latter study to time-dependent injection 
and calculate, by first time, the neutrino output of a hot magnetized corona 
around a galactic black hole. 



2 Static corona model 



2.1 Basic scenario 



The model considered here represents a static corona with a component of 
relativistic particles that can be removed by diffusion (see Romero et al.l 2010l 
for details). In the Bohm regime, the diffusion coefficient is D{E) = Vgc/S, 
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where = E/{eB) is the gyro-radius of the particles. The diffusion rate is 



Relativistic particles injected in the corona can have a local origin. Strong 
shocks resulting from fast magnetic reconnection events and converging mag- 
netic mirrors can in principle accelerate particles up to relativistic energies 

throu gh first-order Fermi mechanism (e.g.. lTsuneta fc Naitolll998l : Ide Gouveia Dal Pino et al, 
20101 ). More complicated processes resulting in second-order Fermi acceler- 



ation are plausible in a turbulent corona (e.g., iDermer et al.l 119961 ). In the 
present work we assume first-order Fermi acceleration. In such the 
acceleration rate for a particle of energy E in a magnetic field B is given by 



rjecB 



(2) 



where < 1 is a parameter that characterizes the efficiency of the acceleration. 
We fix 7] = 10~^, which describes the efficient acceleration by shocks with 
Vs ~ O.lc in the Bohm regime. 

We consider a two-temperature corona in steady state, with a thermal emis- 
sion characterized by a power-law with an exponential cutoff at high energies. 



as ob served in several X-ray binaries in the low-hard state (e.g.. lRomero et al. 



20021 ). Since the corona is in steady state, the assumption of equipartition be- 
tween the different component allows to estimate the mean value of the main 
parameters in the model. Table ([1]) summarizes these values. 

In a corona characterized by such parameters, it is expected that relativis- 
tic electrons and muons lose energy mainly because of synchrotron radiation 
and inverse Compton scattering, whereas for protons and charged pions the 
relevant cooling processes are synchrotron radiation, inellastic proton-proton 
collisions, and photomeson production. Figure ([1]) shows the cooling time for 
the different radiative processes under the conditions of Table ([1]). 



2.2 Spectral energy distributions 



We study the effect of the injection of a power-law distribution of relativistic 
electrons and protons. The main products of hadronic interactions are charged 
pions, which quickly decay producing muons and neutrinos. Neutral pions yield 
gamma-rays, that are a source of secondary pairs. Therefore, we also include 
the effect of all these secondary particles in our treatment. 
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Table 1 

Model parameters 



Parameter 


Value 


Mbh: black hole mass [Mq\ 




Re- corona radius [cm] 


5.2 X W^^'^> 


Te'. electron temperature [K] 




Tj: ion temperature [K] 




Ec- X-ray spectrum cutoff [keV] 


150 


a: X-ray spectrum power-law index 


1.6 


r]: acceleration efficiency 


10^2 


Be- magnetic field [G] 


5.7 X 10^ 


ni,nf.: plasma density [cm~^] 


6.2 X 10^3 


o: hadron-to-lepton energy ratio 


100 


kT: disk characteristic temperature [keV] 


0.1 



Typical value for Cygnus X-1 in the low-hard state 
(2) 35i?G, Rg = 



We solve the transport equation in steady state obtaining particle distribu- 
tions for the different species. Then, the spectral energy distributions (SEDs) 
of all radiative processes are estimated. The results are shown in Fig. ([2]). 

Gamma-rays produced in the corona can be absorbed by different mecha- 
nisms. The most relevant one is photon-photon annihilation. The absorption 
can be quantified by the absorption coefficient or opacity r. In Fig. ([3]) we 
show the opacity due to the interaction between gamma-rays and thermal X- 
ray photons from the corona, which are by far the dominant electromagnetic 
component at low energies. 

Because of the high values of the opacity, it is expected a large number 
of secondary pairs. The effects of internal absorption and the radiation emit- 
ted by secondary pairs are also include in the final SED, which is shown 
in Fig. (jl]). This figure also shows the spectrum of the well-k nown source 
Cygnus X-1, detected by COM PTEL flMcConnell et all l2000h and INTE- 



GRAL (jCadoUe Bel et al.l 120061 ). and the radio emission from the jet. As it 
can be seen, the emission of the corona at low energy is negligible compared 
to that of the jet, in accordance with the idea that both components are 
present in the low-hard state. 
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Fig. 1. Radiative losses in a corona characterized by the parameters of Table [T] 
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Fig. 2. Non-thermal contributions to the total luminosity. 
2.3 Mass Scaling 



Currently, the best evidence for the existence of a common mechanism 
operating in both X-ray binaries and Active Galactic Nuclei (AGNs) is given 
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Fig. 3. Internal absorption due to photon-photon pair production. 
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Fig. 4. Spectral energy distribution obtained with the steady state model of a static 
corona. The pre diction fits the observa t ion made by COMPTE L and INTEGRAL 
of Cygnus X-1 (|McConnell et ahlboool : ICadolle Bel et al.1 12OO6I ') . The non-thermal 



radio emission from the jet is also shown ([Stirling et al 



2001 



by the detectio n of steady jets in accretion regimes with low rates on all scales 
( lMarkofj|2005[ ). It is thought that most black holes spend a significant amount 
of time in this low- luminosity regime on their way in and out of the quiescent 
ground state. For galactic black holes this corresponds to the hard state. In 
AG Ns this w ould correspond to the class of low-luminosity AGNs (LLAGNs; 
e.g.. lHoll2005l) which includes most nearby AGNs such as M81, NGC 4258 and 
M31. 



The latest results strongly support that accretion at low rates is similar 
across the entire range of black hole masses. It seems that the same physical 
model, in which all parameters were expressed in mass-scaling units (such as 
Eddington luminosity and rg), should be able to describe the spectra of either 
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a hard state XRB or a LLAGN, with similar internal parameters. 



In iMarkofj fl2010h is shown that the same outflow-dominated (including 



an outflowing corona as the jet base) model can provide statistically good 
descriptions of broadband data from hard state XRBs as well as LLAGNs, with 
similar ranges in free parameters. These results agree with the idea that the 
physics of at least weakly accreting black holes scales predictably with mass. 
The data are consistent, across the mass scale, with a compact outflowing 
corona directly feeding a continuous, steady, weakly accelerated jet. 

All this suggests that the results presented on this paper can be scaled to 
AGN coronae as well. 



3 Flare model 



3.1 Particle injection 



The temporal dependence of the assumed particle injection is characterized 
by a FRED (Fast Rise and Exponential Decay) behavior, whereas the energy 
dependence is a power-law. The transient injection can be represented by 

Q{E,t) = Q^E-'^e-'^'^-^^^il - e*/^"-) 



^ ^ ''"plat 

— — arctan ( 



(3) 



where Trise = 30 min, Tdcc = 1 h and Tpiat = 2 h. Figure ([5]) shows the primary 
electron injection at = 10'' eV. A similar injection takes place for protons. 
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Fig. 5. Electron injection as a time function at Eq = 10^ eV. 

The power-law has the standard index of a = 2.2. The normalization constant 
Qq can be obtained from the total power injected in relativistic protons and 
electrons, Lrei = Lp + L^. This power is assumed to be a fraction of the 
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luminosity of the corona, Lrci = QrciLc- In the steady state the best fit to 
the observations is obtained with grei = 0.2. During the flare the number of 
relativistic particles increases. In our model, the power injected in the flare 
doubles that of the steady st ate. In so lar flares however, the non-thermal 
injection can be much higher (ILinl l2008l ) suggesting that stronger flares can 
occur in XRB coronae due to massive magnetic reconnection. 



The way in which energy is divided between hadrons and leptons is un- 
known, but different scenarios can be taken into account by setting Lp = aLg. 
We consider a model with a = 100 (proton-dominated scenario, as for galactic 
cosmic rays). 



3.2 Transport equation 



The particle dist ributions N(E,t) can be deriv ed from the solution to the 
transport equation (IGinzburg fc Syrovatskiilll964l ): 



dN{E,t) d 



dt 



+ —{b{E)N{E,t)) + 



N{E,t) 



QiE,t), 



(4) 



where b(E) = ^ . The corresponding solution is: 

loss 



N{E,t) = J Q{E',t- r)e-"(^'^')/*-^dE', 



(5) 



where E^s is implicity defined by 

EeB 



t 



dE' 



(6) 




The primary electron distribution is shown in Fig. ([6]). Similar distributions 
where calculated for all types of particles: p, tt^, 7r°, /x^ and secondary e^. 
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Fig. 6. Evolution in time of electron distribution. 
3.3 Spectral energy distribution 



Figure ([7]) shows the evolution of the total non-thermal luminosity absorbed 
by photon-photon pair production at different energies. The emission at E = 1 
MeV is dominated by the thermal contribution from the corona. 





{h)E 



Fig. 7. SEDs evolution at different energies. The sensitivity of CTA at ii^ = 1 TeV 
IS 2.62 X 10"^'^ TeV cm s , which corresponds to a luminosity of 7.8 x 10 erg 
s~^ at a distance of 2 kpc. Thus the emission might be detectable. 



3.4 Neutrino emission 



The secondary particles produced by hadronic and photohadronic interac- 
tions decay as: 



(8) 
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(9) 



We consider neutrino production by these two c hannels. Thus, the total emis- 
sivity of neutrinos is (IReynoso fc Romeroll2009l ): 



Q,{E, t) = Q^^uiE, t) + Q^^,{E, t), 



(10) 



where 



Qtt^u{E, t) 



dEj-\^^{E^)N^{E^,t) 



e(l-r. 



x] 



EJl-r^) 



(11) 



with X = E/Et^, and 



/dE 
-^t-X{E,)N,^{E,,t) 



-X 



(12) 



In this latter expression, x = E/ E^, /^{i,2} = lA. '^^ ^^^^1 /i{3,4} = /i^ 

The total cosmic ray energy density in a source at a distance d to display 
detectable neutrino emission is Lp/And'^c. Assuming the distance to Cygnus 
X-1 [d ^ 2 kpc), we obtain a cosmic ray energy density of 2.24 x 10~^ eV 
cm~^ at the peak of the emission. 

The differential flux of neutrinos arriving at the Earth can be obtained as: 
1 



,{-.+} 



dE 47rd2 



d'rQ,{E,t). 



(13) 



V 



This quantity, weighted by the squared energy, is shown in Fig. ([8]). Also 
the IceCube sensitivity for 1 year of operation is shown in the figure. Assuming 
that the duty cycle of flares in galactic black holes is around 10 %, the Ice 
Cube detector will be able of detecting neutrinos from a source at ~ 2 kpc 
after 10 years of observations. 



4 Summary and conclusions 



According to our results, the accumulated signal due to neutrino bursts 
in black hole coronae might be detectable for sources within a few kpc. The 
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Fig. 8. Evolution of the neutrino flux, atmospheric neutrino flux and Ice Cube 
sensitivity. 



electromagnetic part of the flare can also be detectable at some energies, de- 
pending of the optical depth for photon annihilation. It should be remarked 
that in some cases the source can be 'dark' at gamm a rays and nonethe less 
a significant neutrino emitter (e.g. Cygnus X-3, see ICerutti et al.l l201ll for 
estimates of the 7 — 7 annihilation) . 



In a future work, we will take into account the effects of absorption due to 
the presence of the anisotropic photon field of the donor star. T his absorption 
depends on the position in t he orbit of the compact object (IRomero et al. 
2OIOI : iRomero fc Vieyrdl201ll ). Another effect that will be studied, is the elec- 
tromagnetic flare modulation produced by the star. This effect is relevant only 
in sources where the duration of the flare is comparable with the period of the 
compact object. Particular emphasis will be put on the source Cygnus X-3, 
which presents a short period of less than 5 hours. 
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